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ABSTRACT 

We present the first combinatorial polynomial time algo- 
rithm for computing the equilibrium of the Arrow-Debreu 
market model with linear utilities. Our algorithm views the 
allocation of money as flows and iteratively improves the 
balanced flow as in [Devanur et al. 2008] for Fisher's model. 
We develop new methods to carefully deal with the flows and 
surpluses during price adjustments. In our algorithm, we 
need 0(n 6 log(n!7)) maximum flow computations, where n 
is the number of persons and U is the maximum integer util- 
ity, and the length of the numbers is at most 0(n\og(nU)) 
to guarantee an exact solution. Previously, [Jain 2007] has 
given a polynomial time algorithm for this problem, which 
is based on solving a convex program using the ellipsoid al- 
gorithm. 

1. INTRODUCTION 

We provide the first combinatorial polynomial algorithm for 
computing the model of economic markets formulated by 
Walras in 1874 12 . In this model, every person has an 
initial distribution of some goods and a utility function of 
all goods. The market clears at a set of prices if every person 
sells their initial goods and then uses their entire revenue to 
buy a bundle of goods with maximum utility. We want to 
find the market equilibrium in which every good is assigned a 
price so that the market clears. In 1954, two Nobel laureates, 
Arrow and Debreu [2], proved that the market equilibrium 
always exists if the utility functions are concave, which is 
why we call this model "Arrow-Debreu market". But, their 
proof is based on Kakutani's fixed point theorem and is non- 
constructive. Since then, many algorithmic results studied 
the linear version of this model, that is, all utility functions 
are linear. 

The first polynomial algorithm for the linear Arrow-Debreu 
model was found by Jain in 2007 [9] ; it is based on solving a 
convex program using the ellipsoid algorithm and simultane- 
ous diophantinc approximation. Before that, Devanur and 
Vazirani [5] gave an approximation scheme for computing 



the Arrow-Debreu model with running time 0(— log — ), im- 
proving [10]. Recently, Ghiyasvand and Orlin [8] improved 
the running time to 0(— (m + nlogn)), where m is the num- 
ber of pairs (i, j) such that buyer i has some utility for pur- 
chasing good j. 

Many combinatorial algorithms consider a simpler model 
proposed by Fisher (see [3]). Eisenberg and Gale [B] reduced 
the problem of computing the Fisher market equilibrium to 
a concave cost maximization problem and thus gave the first 
polynomial algorithm for the Fisher market by ellipsoid algo- 
rithm. The first combinatorial polynomial algorithm for an 
exact linear Fisher market equilibrium is given by Devanur 
et al [4]. They use the maximum flow algorithm as a black 
box in their algorithm. When the input data is integral, their 
algorithm needs 0(n 5 log U + n 4 log e ma x) max-flow compu- 
tations, where n is the number of buyers, U the largest inte- 
ger utility, and e ma x the largest initial amount of money of 
a buyer. If we use the common 0(n 3 ) max-flow algorithm 
(see Q]), their running time is 0(n 8 log U + n 7 log e max ). Re- 
cently, Orlin [11] improved the running time for computing 
the linear Fisher model to 0(n 4 log U + n 3 log e max ) and also 
gave the first strongly polynomial algorithm with running 
time 0(n 4 log n). 



Our results. In this paper, we extend the method in [4] to 
find a combinatorial algorithm for the linear Arrow-Debreu 
market. W.l.o.g., we can assume that each of the n per- 
sons has only one unit of goods, which is different from the 
goods other people have. When the price of that goods is 
Pi, the person will have a budget of p, amount of money 
equivalently. It can be shown that each person only buys 
their favorite goods, that is, the goods with the maximum 
ratio of utility and price. We construct a graph on the nodes 
of buyers and goods as well as equality edges representing 
each buyer's favorite goods, and then find a maximum flow 
in it. In the case of market clearing prices, all the buyers 
and goods will be saturated. There will otherwise be sur- 
plus, so we need to increase the prices of some goods in order 
to decrease the surplus. 

As in [4], we use the notion of "balanced flow", which is a 
maximum flow that balances the surplus of buyers. It is 
helpful for finding a cut in the residual graph. When we 
increase the prices of some goods, the buyers in possession 
of these goods will obtain more budget. Unlike [UJ, where 
the flows stay unchanged during price adjustment, we also 



increase the flows to those goods by the same ratio. So, some 
buyers will have more surplus, and some will have less. We 
can show that in some cases, the new flow is more balanced, 
that is, the /2-norm of the surplus vector will decrease by 
a factor of 1 — f2(l/n 3 ). In other cases, some prices will 
increase by a factor of l + 0(l/n 3 ). Since we can bound the 
number of the latter type of price adjustment iterations by 
a polynomial of n and log U, the total running time is also 
polynomial. 

2. MODEL AND DEFINIITONS 

We make the following assumptions on the model as in Jain's 
paper [5]: 

1. There are n persons in the system. Each person i has 
only one good, which is different from the goods other 
people have. The good person i has is denoted by good 
i. 

2. Each person has only one unit of good. So, if the price 
of good i is pi, person i will obtain pi units of money 
when selling its good. 

3. Each person i has a linear utility function ^ ■ Uij Zij , 
where Zij is the amount of good j consumed by i. 

4. Each Uij is an integer between and U. 

5. For every i, there is a j such that Uij > 0. (Everybody 
likes some goods.) 

6. For every j, there is an i such that Uij > 0. (Every 
good is liked by somebody.) 

7. For every proper subset P of persons, there exist i g P 
and j P such that Uij > 0. 

All these assumptions, with the exception of the last, are 
without loss of generality. The last assumption implies that 
all the equilibrium prices are nonzero [9i , and it is only useful 
for the next section. In Section^ we will discuss more about 
the last assumption. 

Let p = (pi,p2, ■■■,p n ) denote the vector of prices of goods 
1 to n, so they are also the budgets of persons 1 to n. 
In this paper, we denote the set of all buyers to be B = 
{61, 62, b n } and the set of all goods to be C = {ci, C2, c n } 
So, if the price of goods d is pi, buyer bi will have pi amount 
of money. For a subset B' of persons or a subset C of goods, 
we also use p{B') or p{C) to denote the total prices of the 
goods the persons in B' own or the goods in C' . For a vector 
v = {vi,V2, ■••,%), let: 

• \v\ = \vi\ + \v-2\ + ... + \vk\ be the Zi-norm of v. 

• \\v\\ = \/v1 + v\ + ... + vf, be the Z2-norm of v. 

Clearly, every person is happiest with goods j, which max- 
imizes the ratio Uij/pj. Define its bang per buck to be 
an = maxj { /pj}. The classical Arrow-Debreu [2] theorem 
says that there is a non-zero market clearing price vector. 



3. THE ALGORITHM 

Our algorithm resembles [4., which finds a balanced flow 
and increases the prices in the "active subgraph". But, in 
the Arrow-Debreu model, when we increase the prices of 
some good i, the budget of buyer i will also increase. So, we 
need to find a careful way to prevent the total surplus from 
increasing. 

Construct a flow network G = ({s, t}uBL)C, Eg), where s is 
the source node and t is the sink node, then B = {61, b n } 
denotes the set of buyers and C — {ci, ...,c„} denotes the 
set of goods. Eg consists of: 

• Edges from s to every node bi in B with capacity pi . 

• Edges from every node Ci in C to t with capacity pi. 

• Edges from bi to c, with infinite capacity if Uij/pj — 
Qi. Call these edges "equality edges". 

So, our aim is to find a price vector p such that there is 
a flow in which all edges from s and to t are saturated, 
i.e., (s,C U B Ut) and (sUC U B ,t) are both minimum cuts. 
When this is satisfied, all goods are sold and all of the money 
earned by every person is spent. 

In a flow /, define the surplus r(bt) of a buyer i to be 
the residual capacity of the edge (s, 6^), and define the sur- 
plus r(cj) of a good j to be the residual capacity of the 
edge (cj,t). That is, r(6 4 ) = p t - Yljfiii and r i c j) = 
Pj — Y^,i fijj where fij is the amount of flow in the edge 
(bi,Cj). Define the surplus vector of buyers to be r(B) = 
(r(bi), ^(62), r(b n )). Also, define the total surplus to be 
|r(B)| = X^ r (fr*)i which is also ^ . r(cj) since the total 
capacity from s and to t are both equal to F° r con ~ 

venience, we denote the surplus vector of flow /' by r'(B). 
In the network corresponding to market clearing prices, the 
total surplus of a maximum flow is zero. 

3.1 Balanced flow 

As in [I], we define the concept of balanced flow to be a 
maximum flow that balances the surpluses of buyers. (How- 
ever, unlike in their paper, the surpluses of goods can be 
positive here, which are not supposed to be balanced, so the 
balanced flow is not necessarily unique.) 

Definition 1. In the network G of current p, a balanced 
flow is a maximum flow that minimizes ||r(B)|| over all 
choices of maximum flows. 

For flows / and /' and their surplus vectors r(B) and r'(B), 
respectively, if ||r(B)|| < ||r'(B)|], then we say / is more 
balanced than /'. The next lemma shows why it is called 
"balanced". 



Lemma 1. /^/ If a > bi > 0,i = 1,2, ...,k and 5 > 
X^ - =i where 8, Si > 0, i = 1, 2, k, then: 

||(a,6i,6 2 ,...,fefc)|| 2 < \\a + 5,b 1 -5 1 ,b 2 -5 2 , ...,b k -5 k \\ 2 -5 2 . 



Proof. 
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they disappear from the equality graph. If there are edges 
(h,Cj) and 0,c fe ) where b z G B(S),Cj,c k G T(B(S)), then 
U ij/Pj = Uik/pk- Since the prices in T(B(S)) are multiplied 
by a common factor x, Uij/pj and Uik/pk remain equal after 
a price adjustment. However, the goods in C\T(B(S)) will 
become more attractive, so there may be edges from B(S) to 
C\F(B(S)) entering the network, and the increase of prices 
needs to stop when this happens. Define such a factor to be 
X(S), that is, 

X(S) = min{^-^\b t G B(S), (b t ,c 3 ) G E G ,c k £ T{B(S))}. 

Pj Uik 



Lemma 2. Jjjj In the network G for a price vector p, given 
a maximum flow f , a balanced flow /' can be computed by 
at most n max-flow computations. 

PROOF. In the residual graph Gf w.r.t. to /, let S C 
B U C be the set of nodes reachable from s, and let T — 
(B U C) \ S be the remaining nodes. Then, there are no 
edges from S n B to T n C in the equality graph, and there 
is no flow from TnBto5nC. The buyers in T n B and the 
goods in S H C have no surplus w.r.t. /, and this holds true 
for every maximum flow. Let G' be the network spanned 
by s U S U t, and let /' be the balanced maximum flow in 
G' . The /' can be computed by n max-flow computations. 
(Corollary 8.8 in [¥] is applicable since (s U S, t) is a min- 
cut in G' .) Finally, /' together with the restriction of / to 
s U T U t is a balanced flow in G. □ 



The surpluses of all goods in /' are the same as those in / 
since we only balance the surplus of buyers. 

3.2 Price adjustment 

As in [H [TT] , we need to increase the prices of some goods to 
get more equality edges. For a subset of buyers Bi, define 
its neighborhood F(_Bi) in the current network to be: 

r(-Bi) = {Cj G C\3b t G Bi,s.t.(bi,cj) G E G }- 

Clearly, there is no edge in G from B\ to C \ T(B\). In 
a balanced flow /, given a surplus bound 5* > 0, let B(S) 
denote the subset of buyers with surplus at least S, that is, 
B(S) = {bt G B\r{bi) > S}. 



Lemma 3. In a balanced flow f, given a surplus bound S, 
there is no edge that carries flow from B\B(S) to T(B(S)). 

Proof. Suppose there is such an edge (bi,Cj) that car- 
ries flow such that h £ B(S) and c, G F(B(S)). Then, 
in the residual graph, there are directed edges (bk,Cj) and 
(cj,bi) with nonzero capacities in which b k G B(S). How- 
ever, r(bk) > S > r(bi), so we can augment along this path 
and get a more balanced flow, contradicting that / is already 
a balanced flow. □ 



From Lemma[3] we can increase the prices in T(B(S)) by the 
same factor x without inconsistency. There is no edge from 
B(S) to C\T(B(S)), and the edges hom B\B{S) to T(B(S)) 
are not carrying flow, and hence, there will be no harm if 



So, we need 0(n 2 ) multiplications/divisions to compute X(S). 
When we increase the prices of the goods in T(B(S)) by a 
common factor x < X(S), the equality edges in B(S) U 
T(B(S)) will remain in the network. We will also need the 
following theorem to prevent the total surplus from increas- 
ing. 

Theorem I. Given a balanced flow f in the current net- 
work G and a surplus bound S, we can multiply the prices 
of goods in F(B(S)) with a parameter x > 1. When x < 
xmm{ Pi /(s)i - r(bi))\bi G B{S),a £ V(B(S))} and x < 
X(S), we obtain a flow /' in the new network G' of adjusted 
prices with the same value of total surplus by: 



f' 



x-fa ifc 3 eT(B(S)); 

fa iS * tv(B{S)). 



Then, the surplus of each good remains unchanged, and the 
surpluses of the buyers become: 

x-r(h) ifbi€B{S),Ci€T(B(S)); 

(1 - x) Pi + x ■ r(bi) if h G B(S), a i T{B{S)); 

(x - l) Pi + r(bi) if h £ B(S), a G T(B(S)); 

r(bi) ifb l £B{S),c l £T{B{S)). 

We call these kinds of buyers type 1 to type 4 buyers, respec- 
tively. 

Proof. Since the flows on all edges associated with goods 
in F(B(S)) are multiplied by x, the surplus of each good in 
r(B(S)) remains zero. Only the surplus of type 2 buyers 
decreases because the flows from a type 2 buyer bi are mul- 
tiplied by x, but its budget pi is not changed. The flow after 
adjustment is x(pi — r(bi)). We need this to be at most pi, 
so x < Pi/(pi — r(bi)) for all type 2 buyers bi, and in /', the 
new surplus r'(6<) = (1 — x)pi + xr(pi). 

Since both money and flows are multiplied by x for a type I 
buyer, his surplus is also multiplied by x. For a type 3 buyer 
bi, his flows are not changed, but his money is multiplied by 
x, so the new surplus is xpi — (p» — r{bi)). □ 

After each price adjustment, in the new network, we will 
find a maximum flow by augmentation on the adjusted flow 
/' and then find a balanced flow by Lemma This will 
guarantee that when the surplus of a good becomes zero, it 
will not change to non-zero anymore. Thus, the prices of 
the goods with non-zero surplus will not be adjusted. 

Property 1. The prices of goods with non-zero surpluses 
remain unchanged in the algorithm. 



3.3 Whole procedure 

The whole algorithm is shown in Figure [1] where A" is a 
constant we will set later. In this section, one iteration de- 
notes the execution of one entire iteration inside the loop. 
We will discuss the rounding and termination conditions in 
Section EH 



In the first iteration, we constructed a balanced flow / in 
the network where all prices are equal to 1. In the equality 
graph, we have at least one edge incident to every buyer. 
The total surplus will be bounded by n, actually n — I as at 
least one good will be sold completely. From Theorem [T] in 
the execution of the algorithm, the total surplus will never 
increase. 

To ensure that the algorithm will terminate in a polynomial 
number of steps, we will require the following lemmas. From 
Property [T] the prices of goods with surplus stay one during 
the whole algorithm, so there is still a good with price one 
in the end. And, we need to bound the largest price: 



Lemma 4. The prices of goods are at most (nU) n . 

Proof. It is enough to show that during the entire al- 
gorithm, for any non-empty and proper subset C of goods, 
there are goods Cj £ C,Cj ^ C such that Pi/pj < nil. So, 
when we sort all the prices in decreasing order, the ratio of 
two adjacent prices is at most nil. Since there is always a 
good with price 1, the largest price is < (n[/) n_1 . 



• be multiplied by a factor of 1 + 0(l/n 3 ) when x — 
1 + -Js, or 

• be divided by a factor of 1 + f2(l/n 3 ). 

Note that by Lemma [S] the rounding procedure can only 
increase ||r(_B)|] by a factor of l + 0(l/n 4 ) since S > e/(e-n). 
(We will leave the discussion of Lemma [5] later, so we can 
ignore it in the analysis.) 



Theorem 2. In total, we need to compute 0(n 6 log(n(7)) 
maximum flows, and the length of numbers is bounded by 
0(n log(n(/)) . Thus, if we use the common 0(n 3 ) max-flow 
algorithm (see \Tj), the total running time is 0(n 10 log 2 (n(7)). 

Proof. By Lemma f3] every price can be multiplied by 
x = 1+ ^ for OQog 1+1/Kn3 (nU) n ) = 0(n 4 log(n[/)) 
times, so the total number of iterations of the first type 
is 0(n 5 log(n[/)). The total factor multiplied to ||r(S)|| by 
the first type iterations is (1 + 0(l/n 3 )) 0( ™ 5 log(nC/)) . 

At the beginning, |]r(_B)|| < ~Jn. When the algorithm ter- 
minates, ||r(_B)|| < e — An i v i n , so the number of second 
type iterations is bounded by 

log 1+ n(i Al 3 ) (V(l + 0(l/n 3 ))°(" 5 W)) 
= 0(n 5 log(nf/)). 



If C contains goods with surpluses, then their price is 1. The 
claim follows. 

Let B — r(C) be the set of buyers adjacent to goods in C 
in the equality graph. If there exist bi, Cj s.t. bi £ B, Cj ^ C 
and Uij > 0, let Ck £ C be one of the goods adjacent to 
bi in the equality graph, and then Uij/pj < Uik/pk- So, 
Pk/Pj < v-ik/Uij < U. 

If there do not exist such bi,Cj, then there is no flow between 
B and C \ C, and there is bk ^ B, but c& £ C. Otherwise 
the persons whose goods are in C will not like any goods 
not in C, contradicting assumption (7). Let B' — {j\bj £ 
B, Cj £ 6} and B" = {j\b 3 £ B, cy £ C}. We have: 

Pk < p(B") = p(C) - pdjlbj £ B, Cj £ C}) 

< p(B) - P ({j\b 3 eB, C] eC}) 

= p(B'). 

The inequality of the second line holds since goods in C 
have surplus and all of the flows from B go to C. Thus, 
B' must be non-empty, and hence, there is a j £ B' with 
Pi > p(B')/n. We conclude pk < npj. □ 

By Lemma [9] we can round to the exact solution when the 
algorithm terminates. To analyze the correctness and run- 
ning time, we need the following lemma: 

Lemma 5. After every price adjustment by x, the li -norm 
of the surplus vector ||r(B)|| will either 



Thus, the total number of iterations performed is bounded 
by 0(n 5 \og(nlI)). Since we need to compute n max-flows for 
the balanced flow in every iteration, we need 0(n 6 log(n!7)) 
maximum flow computations in total. By Lemma [4] and 
Lemma [S] the prices are rational numbers < (nt/) n_1 and 
with denominator < U n A = in 9 U 4n . Thus, the length of 
the numbers to be handled is bounded by 0(n \og(nU)). 
Note that max-flow computations only need additions and 
subtractions. We perform multiplications and divisions when 
we scale prices and when we set up the max-flow compu- 
tation in the computation of balanced flow. The numbers 
of multiplications/divisons is by a factor n less than the 
numbers of additions/subtractions, and hence, it suffices to 
charge 0(nlog(nU)) per arithmetic operation. □ 



Next we will prove Lemma [5] When we sort all the buy- 
ers by their surpluses 61,62,— ,b„, 61 is at least \r(B)\/n 
(where |r(i?)| is the total surplus). So, for the first i in which 

rl^ll) > 1 + V", we can see < 1 + 1/n for j < i, 

so r(6i) > r(6i)(l + l/n)~ n > |r(B)|/(e ■ n). When such an 
i does not exist, each r(bi) is larger than \r(B)\/(e ■ n), and 
all goods in T(B) must have zero surplus because the flow 
is otherwise not maximum. Thus, there are goods that have 
no buyers, and hence, either new equality edges emerge, or 
x reaches 1 + ■ 



(condition (3a) below). 



From the algorithm, in every iteration, x satisfies the follow- 
ing conditions: 



Initially set p; = 1 for all goods i\ 
Repeat 

Construct the network G for the current p, and compute the balanced flow / in it; 
Sort all buyers by their surpluses in decreasing order: bi, 62, b n ~, 
Find the first i in which r ^ ft '^ > 1 + 1/n, and i = n when there is no such i; 
Let the surplus bound S = r(6») and obtain B(S),F(B(S)),X(S); (B(S) = {bi, b 2 , bi}) 
Multiply the prices in T(B(S)) by a gradually increasing factor x > 1 until: 
(Let /' be the flow corresponding to x which is constructed according to Theorem [T]) 
New equality edges emerge (x reaches X(S)); 

OR the surplus of a buyer £ B(S) and a buyer ^ B(S) equals in /'; 
OR x reaches 1 + „ u 
Round the prices in F(B(S)) according to Lemma [8] with A = 4n 9 U ; 
Until \r(B)\ < e, where e = An i v $ n ; 

Finally, round the prices according to Lemma [U] to get an exact solution. 



Figure 1: The whole algorithm 



(2) In /', r'(6) > r'(b') for all b £ B(S),b' $ B{S). Here, 
r'(b) is the surplus of b w.r.t. /', the flow corresponding 
to x by Theorem Q] 

(3) If x < 1 + ^rr, the following possibilities arise: 

(a) There is a new equality edge (bi,Cj) with 6^ g 
B(S),Cj g F(B(S)). By Lemma |B] below, we can 
obtain a flow /" in which either r" (bi) = r'(bi)—pj, 
or there is a bk B(S) with r"(bi) = r"(bk) (same 
as (b)). 

(b) When x satisfies the second requirement in the al- 
gorithm, it satisfies: there exists b G B(S) and 
b' £ B(S) such that r'(b) = r'(b') in /'. 



unchanged for type 4 buyers. Note that the surplus of a 
type 1 or 2 buyer cannot be smaller than the surplus of any 
type 3 or 4 buyer. From Theorem [T] and condition (2), we 
infer that the total surplus will not increase, type 2 and 3 
buyers will get more balanced, and r'(b) — x ■ r(b) for type 
1 buyers b, so \\r'(B)\\ < x\\r(B)\\ = (1 + 0(l/n 3 ))\\r(B)\\. 



In (3a), there is a new equality edge (bi 



After the pro- 



cedure described in Lemma [6] if there is no bk ^ B(S) such 
that r"(h) = r"(b k ), then r"(bi) = r' (h) - Pj ( Pj > 1). For 
all b k £ B(S), r"(bi) > r"(bk), and r"(b k ) = r'(b k ) + S k , 
where 8 k > and X^^sfs) ^fc — Pj- Because |r(B)j < n, 
\\r(B)\\ 2 < n 2 . By Lemma|T] 



Lemma 6. If there is a new equality edge (bi,Cj) with bi £ 
B(S),Cj (f: T(B(S)), we can obtain a flow f" from /' in 
which either r" ' (bi) = r'(bi) — Pj, or there is a b k B(S) 
withr"(bi) = r"(b k ). 

Proof. Let B' C B \ B(S) be the set of buyers with 
flows to Cj in /', and let w be the largest surplus of a buyer 
in B \ B(S). Run the following procedure (/" denotes the 
current flow in the algorithm): 



!k"(B)|| 2 < ||r'(B)|| 2 - Pj 2 

< a; 2 ||r(B)|| 2 - 1 

< x 3 ||r(B)|| a --^||r(B)|| : 
= (l-e(l/n 2 ))||r(B)|| 2 . 



Augment along (bi,Cj) gradually until: 

r"(h) =i»or r"( Cj ) = 0; 
If r"(bi) = w then Exit; 
For all b k G B' in any order 

Augment along (bi,Cj,bk) gradually until: 
r"(bi) = max{r"(fe fe ),™} or /"(&„, cj) = 0; 

Set w = max{r"(6fe), w}; 

If r"(bi) = w then Exit. 

During the procedure, the surplus of bi decreases but cannot 
become less than the surplus of a buyer in B \ B(S), so 
condition (2) holds. In the end, if r"(bi) = w, then there 
is a bk 6 B \ B(S) s.t. r"(bi) = r"(bk); otherwise, Cj has 
no surplus, and the flow to it all comes from bi, so r" (bi) = 
r'(bi)~ Pj . □ 

From Theorem [TJ the surpluses in /' will increase for type 1 
and 3 buyers, will decrease for type 2 buyers, and will stay 



So, we have ||r"(B)|| = (1 - fl(l/n 2 ))\\r(B)\\. 

In (3a), after the procedure described in (3a), if there exists 
bk ^ B(S) such that r" (bi) = r"(bk), then we are in a similar 
situation as (3b), possibly with an even smaller total surplus. 
So, we can prove this case by the proof of (3b). 

In (3b), let ui, 112, ut and vi,V2, vy be the list of orig- 
inal surpluses of type 2 and 3 buyers, respectively. De- 
fine u = min{ui},v = max{vj}, so Ui > u for all i, and 
Vj < v for all j, and u > (1 + l/n)v. After the price and 
flow adjustments in Theorem [T] the list of surpluses will be 
mi — Si,u 2 — S 2 , uk — S k and vi + S' ly v 2 + 5' 2 , v k + 5' k , 
(here 5i,8'j > for all i,j), and there exist I, J such that 
ui — Si = vj + 5'j, where ui~5j is the smallest among Ui — 5i, 
and vj + S'j is the largest among Vj + Sj by condition (2). 
Since the surpluses of type 1 edges also increase, we have 
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< —(u — v) 5j 

i 

< — (it — t>) max{<5/, 5j} 



2(n + l) 2 



Lemma 7. In a connected component \& containing k goods 
in the equality graph, if we know that pj is a rational number 
with denominator N, where Cj (E \f r D C, then all the prices 
of goods m $ n C are rational numbers with denominator 
< N ■ U k . 



Proof. Find a tree that connects all the goods in 
The tree will contain k + k' — 1 edges if it contains k! buy- 
ers. Then, we can get k — 1 linear independent equations 
Pj/uij = Pji/uiji when both (6i,c 3 ) and (&j,c,-/) are tree 
edges. Together with the equation pj = 7/JV for some in- 
teger I, we can see that all the prices of goods in have 
denominator < N ■ U k . □ 



Given an integer A, we call a connected component in the 
equality graph consistent if it has a good whose price is a 
rational number with denominator A. Then, by Lemma [7] 
the prices of goods in a consistent connected component are 
rational numbers with denominator < A • U n . 



Let wi,W2, ...wyi be the list of surpluses of type 1 buyers; 
all of them are < e • u. After price adjustment, the surpluses 
will be x ■ Wi,x ■ W2, —x ■ Wk" from Theorem [T] Compute: 



( 1 + ^) 2 £" 



^ E-? + (^3 + TO)-- eV 

i 

\ ^ 2/2 1 \ 2 2 

= > Wi + ( + r je w . 



Let K = 32e 2 , then the change to the sum of squares of 
surpluses for type 2 and 3 buyers is less than — -^u 2 = 
— -^ I e 2 u 2 . The total change to ||r(_B)|| 2 is: 



2 1 

+ 
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\ 2 2 

)e u . 



Since u > ^r(&i) for all buyers &i, nu 2 > -^i\\r{B)\\ 2 . Since 
the change is negative, we have: 



-'{B)\\ 2 < ||r(B)|| 2 + (- 

= lk(B)ll 2 - 



+ ^k)-\Hb)\ 



Kn 2 K 2 n 5 n 



Kn 3 ' 
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Thus, Lemma [5] is proved. 



3.4 Rounding and termination condition 

In this section, we will show how to round the prices to ra- 
tional numbers with denominators of length 0(n log(n[/)). 
Also, we need the rounding process to obtain an exact mar- 
ket equilibrium when the surplus is very small. Here, we 
define the equality graph F on B U C of undirected equal- 
ity edges between buyers and goods, and we consider every 
connected component in this undirected equality graph. 



Lemma 8. In a balanced flow f in the network, given 
A > nU n and a surplus bound S > n /A, if all the con- 
nected components in (B\B(S),C\T(B(S))) are consistent, 
we can adjust the prices in T(B(S)) so that all connected 
components in the equality graph are consistent. In the ad- 
justed flow f by Theorem^ \\r'(B)\\ = \\r(B)\\ (l + 0(^)), 
where r'(B) is the surplus vector in /'. 



Proof. The procedure is shown below: 



Set B' = B(S); 
Repeat 

Multiply the prices in P(-B') by x > 1 until: 
A price in P(-B') has denominator A; 
OR a new equality edge emerges; 
Update the flow / by Theorem [T] 
Remove new consistent components from B' U T(B') 
Until B' = 0. 



Since all the prices change by at most 1/A, the changes to 
the total surplus of type 2 buyers is at most n/A < S, so 
the surplus of a type 2 buyer is still positive. In addition, 
the surplus added to each type 1 or 3 buyer is at most 1/A, 
so 

\\r'(B(S))\\ 2 

< \\r(B(S))\\ 2 + ^\r(B(S))\ + ^ 

< \\r(B(S))f + JL|| r (B(S))|| a + 
= ||r(B(S)||(l + 0(-L)). 



" \HB(S))f 



During the algorithm, we can see that all the connected com- 
ponents in (B \ B' , C \ T(B')) are consistent since we move 
the new consistent components to it. When we find new 
equality edges connecting B' and C \ T(B'), some nodes in 
B'UT(B') will connect to (B\B' ,C\T(B')), so these nodes 
can be removed. When a price in r(-B') has denominator 



A, the component containing it will become consistent, so 
the loop will run for at most n times. In each loop, we need 
to compute 0(n 2 ) multiplications/divisions, so the running 
time for this rounding procedure is less than the computa- 
tion of a balanced flow. □ 



Lemma 9. When the total surplus is < 



a flow f , we can obtain an exact solution from the current 
equality graph. 



Proof. Add the edge (bi,Ci) for each person i to the 
equality graph F to obtain F' . For a connected component 
of F', the sum of the prices on both sides are the same. 
For every component <E> of F' with no surplus node, increase 
its prices by a common factor until a new equality edge 
emerges; this will unite two components. Repeat this until 
all components in F' have a surplus node. We may assume 
w.l.o.g. that F' becomes connected by this process. Other- 
wise, the following argument can be applied independently 
to each component of F' . The total surplus is still less than 
e. The following rounding procedure will be performed on 
these revised prices. 

Denote the set of connected components in the equality 
graph (after the adjustments of the previous paragraph) by 
A = {'I'fe}. For each component ^ in F, find a spanning 
tree T k in it, then write the following equations: 

Pj/uij =Pj'/uijr,V(bi,Cj), (bijCji) G T fe . 

Since we have one such equation if Cj and cy are connected 
by one bi, we can have | ^ n Cj — 1 linear independent equa- 
tions. The total number of linear independent equations for 
all components in F is n — |A|. 

Since there is no flow between components, for each com- 
ponent $fe in F, the money difference between buyers and 
goods in 'l/fc is only the surplus difference. So, we can write 



Consider the following n linear equations with ek removed: 

Pj l u ij = Pj> I u ij' , V (bi,Cj), (h , cy) £ T fc 
J2 E ^ = ' V * fe 



p'i = l, 3 r(a) > 0. 

They can be denoted by Ap' = X' , so there is also a unique 
solution. The solution will be rational numbers with a com- 
mon denominator D < nU n by Cramer's rule. Since \ \X\ — 
\X'\\ < 2e, the difference \p'i — Pi\ of solutions of each price is 
at most 2e • nll n = 2n i u ^ n by Cramer's rule. The difference 
between any two numbers of denominators D, D' < nil 71 is 
a rational number of denominator D • D' < n 2 U 2n , which 
is larger than 2\p\ — Pi\. Since p\ is a rational number with 
denominator D < nU n , we can get p'i by rounding p; to 
the nearest rational number of denominator < nU n . This 
can be done by continued fraction expansion, which needs 
0(nlog 2 D) = 0(n 3 log 2 (nil)) time by Theorem 3.13 in [7J. 
We can also compute D = det(A) directly and round every 
price to the nearest rational with denominator D or solve 
the linear equations Ap' = X' . By Theorem 5.12 in [7], 
computing the determinant of a matrix of dimension n with 
entries < U takes 0(n 4 log U) time, and solving Ap' — X' 
also takes 0(n 4 log U) time. 

Now, all the prices p'i are of the form qi/D, where qi, D are 
integers and D < nil 71 a common denominator. So, \pi — 



= 4=;, in which e — 



D 



Con- 



struct the flow network G' for the new prices q = [qi, q2, q n 
Consider any bi G B and c-j,Ck G C and assume Uij/pj < 
Uik/Pk- Then, 

Uijq k < Uij(p k D + e) 

< UikPjD + Uije' 

< Uik(qj + e ) + ttye 

< u ik qj + {uik + Uij)e 

< Uikqj + 1, 

and hence, Uijqt < uiuqj since Uijqk and Uikqj are integral. 
We conclude that the edges in G are all in G' . 



^2 Pi- ^2 P* = e fc>V* fc . 

Here, et (positive or negative) comes from the surpluses of 
goods and buyers, so ^ |efc| < 2e. If bi and a belong to 
distinct connected components tyj and fyk, the coefficient 
of pi is +1 in the equation of ^j, —1 in the equation for 
V&fe, and in all other equations. If bi and d belong to the 
same connected component, the coefficient of pi is zero in all 
equations. Assume now that there is a proper subset of the 
equations that is linear dependent. Then, if bi or a belongs 
to one of the components in the subset, both of them do. 
However, the subset of components is a proper subgraph of 
F' , and hence, there is at least one i such that only one of 
bi or a belongs to the subset of components. Thus, we have 
|A| — 1 independent equations. 

Since there is a good a with non-zero surplus, we have pi = 
1. Thus, the current price vector p is the solution of these 
linear equations Ap — X in which A is invertible. 



Denote the size of the cuts (s, B U C U t) and (s U B U C, t) 
in G' by Z, which is an integer. Then, the size of this cut 
in G is > (Z — ne)/D. If there is another cut in G' of 
size < Z — 1, it is also a cut in G, and its size in G is 



Z/D 



-2ne' 



D 
l-3e' 



so the maximum flow in G 



will have total surplus > 1 ~^ e > e. Thus, (s,BuCUi) and 
(sUBUC,i) are both min-cuts in G' , so the prices reach a 
market equilibrium. □ 

4. GENERAL CASE 

Here, we consider the case which does not satisfy the as- 
sumption (7) in Section [2] We draw the linking graph of 
persons in which there is a directed edge from i to j iff 
Uij > 0. If the graph is strongly connected, then the case 
satisfies assumption (7). Otherwise, we can shrink each con- 
nected component into one vertex, then the graph will be a 
DAG, and we can find a topological order of strongly con- 
nected components: Pi, P% % ...Pk, in which there are only 
edges from a lower order to a higher order. We use the 
algorithm in Section [3] to compute the equilibrium for all 



the persons in every strongly connected component Pi (i — 
1,2, ...,k). For i — 2,...,k, multiplying the prices in Pi by 
(U + 1) ■ max{pj\j £ Pi-i} will ensure that there are no 
equality edges from Pi to Pj for i < j. Since the persons 
in Pj do not like any goods in P t for i < j, this will not af- 
fect the equilibrium of every component, so we get a global 
equilibrium. 
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